Interaction between impurities and solitons in quasi one-dimensional spin-Peierls 

systems 



D. Augier a , P. Hansen 6 , D. Poilblanc a , J. Riera a ' fc , and E. S0rensen a 
a Laboratoire de Physique Quantique & Unite Mixte de Recherche CNRS 5626 
Universite Paul Sabatier, 31062 Toulouse, France 
b Instituto de Fisica Rosario, Consejo Nacional de Investigaciones Cientificas y Tecnicas y Departamento de Fisica 
Universidad Nacional de Rosario, Avenida Pellegrini 250, 2000-Rosario, Argentina 

(February 7, 2008) 

The role of the spin-phonon coupling in spin-Peierls chains doped with spin-0 or spin-1 impurities 
is investigated by various numerical methods such as exact diagonalization, quantum Monte Carlo 
simulations and Density Matrix Renormalization Group. Various treatments of the lattice, in a 
fully quantum mechanical way, classically in the adiabatic approximation or using a fixed three- 
dimensional dimerization pattern are compared. For an isolated chain, strong bonds form between 
the two spin-| sites next to the impurity site, leading to the appearance of magneto-elastic solitons. 
We also show that these excitations do not bind to spin-0 impurities but are weakly attracted by 
spin-1 impurities. However, the interchain elastic coupling generates an effective confining potential 
at the non-magnetic impurity site which can lead to the formation of soliton-impurity bound states. 
We also predict that a soliton and an antisoliton bound to two impurities on the same chain can 
annihilate each other when the separation between the impurities is smaller than a critical value 
depending on the interchain elastic constant. 

PACS: 75.10 Jm, 75.40.Mg, 75.50.Ee, 75.30.Hx 



I. INTRODUCTION 

Quasi one-dimensional (ID) quantum antifcrromag- 
nets can exhibit surprising magnetic properties at 
low temperature. The recent discovery of the spin- 
Peierls (SP) transition in the inorganic compounds 
CuGe0 3 (Ref. |l|) and NaV 2 5 (Ref. |) has drawn both 
experimental and theoretical interest. The chemistry of 
these compounds indeed allows for the synthesis of large 
single crystals and consequently the achievement of new 
experimental studiesfj which were not accessible to the 
previously known organic SP materials. Nevertheless, 
these two compounds have quite different behaviors pre- 
sumably due to the fact that while NaA^Os is believed to 
have quarter-filled ladders, CuGe03 is well described by 
weakly coupled spin-i Heisenberg chains. In this sense, 
the present studies are more related to this latter com- 
pound. 

At a critical temperature, T c , the SP compounds un- 
dergo a phase transition driven by the spin-phonon cou- 
pling. This SP transition, which is characterized by 
the opening of a spin gap and a lattice dimerization, 
is experimentally signaled by an isotropic drop of the 
magnetic susceptibility revealing the non-magnetic na- 
ture of the ground state (GS). A general physical pic- 
ture can be drawn from the consideration of the exact 
GS of the frustrated Heisenberg chain at the so-called 
Majumdar- Ghosh point. The spontaneously dimerized 
non-magnetic GS is two-fold degenerate, corresponding 
to two possible dimerization patterns A and B, which are 
a succession of disconnected singlet dimers. The elemen- 
tary excitation called soliton (antisoliton) consists of aia 
unpaired spin separating A and B (B and A) patterns!! 



Solitons and antisolitons propagate, then acquiring a dis- 
persion. 

For temperatures above the critical temperature T c , 
quasi- ID SP compounds are usually described as inde- 
pendent uniform AF Heisenberg chains, in some cases 
including terms describing the frustration present in the 
system. The nearest neighbor J and next-nearest neigh- 
bor a J exchange integrals can-then be estimated by a fit 
of the magnetic susceptibilitytffl. The high temperature 
behavior is governed by J and the position of the maxi- 
mum by the frustration ratio a. Values such as J ~ 160 K 
and a « 0.36 have been proposed for CuGeOs (Ref. ||). 

CuGe03 can be easily doped by substituting magnetic 
Ni 2+ or non-magnetic Zn 2+ impurities to spin-i Cu 2+ . 
A rapid suppression of the spin gap due to impurity, dop- 
ing has been measured by magnetic susceptibilitjfl and 
inelastic neutron scatteringa experiments. Competition 
between the antifcgmmagnetic (AF) and SP phases has 
also been observedQ'ErtJ. Theoretical work analyzing the 
effecls-due to the doping by vacancies has been carried 
outlitlij, assuming a static dimerization. 

The coupling to the lattice plays a major role in SP 
compounds. A description of these systems in terms of 
a static dimerized Heisenberg chain is widely found in 
the literature. This approach has some drawbacks: the 
dimerization is introduced de facto in the model and is 
totally frozen, the lattice cannot adjust to spin fluctu- 
ations, which are essential to understand the behavior 
under high magnetic fields, or to inhomogeneities intro- 
duced by impurities. 

In this paper, SP spin-lattice models with spin-0 or 
spin-1 impurities have been studied at T — or very 
low temperature using various numerical techniques such 
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as Exact Diagonalization (ED)0, Quantum Monte-Carlo 
(QMC) and Density Matrix Renormalization Group 
(DMRG). Emphasis has been put on the appearance of 
magneto-elastic excitations and their interplay with im- 
purities. A strictly ID dynamical spin-lattice model is 
first studied in Sec. [□]. Non-uniform dynamical or adia- 
batic lattice distortions are considered and the existence 
of bound states due to impurities is investigated. The 
role of the three-dimensional character of the phonons is 
then considered, first in the simplest description in terms 
of a static dimerized Heisenberg chain. Bound states be- 
tween a soliton and a non-magnetic impurity are quan- 
titatively studied in Sec. III. A simple physical picture 



is al so g iven to explain the spectrum of bound states. In 
Sec. IV a more involved model is introduced by consid- 
ering explicitly the elastic coupling between the ID spin- 
phonon chains. This model enables us to understand the 
physical origin of the soliton-impurity bound states. 



II. ID SPIN-PHONON CHAIN 
A. Models 

The key ingredient of the physics of spin-Peierls com- 
pounds is the dynamical spin-phonon coupling. A spin-i 
excitation is expected to locally distort the lattice, cre- 
ating an elastic soliton. Our starting point is the frus- 
trated Heisenberg chain with the parameters J and a 
defined below. In addition, we assume a linear depen- 
dence of the nearest neighbor exchange integral on the 
relative atomic displacement. Dispcrsionless phonons of 
frequency f2|-are_|Considered for sake of simplicity. One 
then obtainfilTEl, 

< )= JE(i+# ! +t, t ))(s.'4i- \) 

i 

+ ajJ2SrS, + 2 + nJ2b\b l , (1) 

3 i 

where b\ (b i ) is the phononic creation (destruction) oper- 
ator at site i and g is the magneto-elastic coupling. Note 
that we have assumed that the magneto-elastic coupling 
associated with the next-nearest neighbor bonds can be 
neglected. 

The adiabatic limit of Hamiltonian (Q) is of great inter- 
est. When the phonon frequency £1 is sufficiently small 
compared to the spin fluctuation energy scale J, the 
phonon degrees of freedom can be treated classically. The 
Hamiltonian is then written as, 

/(a) 
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where Si is the distortion at site i, K the string constant 
(= £l/(2g 2 )) and the last term corresponds to the elastic 
energy loss. 5i is expressed in units so that the magneto- 
elastic coupling, which can be absorbed in the definition 
of K, is set to unity. It is interesting to notice that the 
frustration a alone (i.e. g = 0) can lead to a dimer- 
ized state for a > a c ~ 0.241167 (Ref. ||(2|g,||,|2|) . 
In fact, the relevance of a is two fold; (i) it is required 
to explain qualitatively the properties of the real com- 
pounds (a ~ 0.36 has been proposed in CuGe0 3 ) and 
(ii) it can be used in numerical simulations to reduce the 
spin-spin correlation length (for a = 0.5 and g = 0, it 
is only one lattice unit) and hence to reduce finite size 
effects. 

The substitution of a spin-i by a spin-S (S=0 or 1) im- 
purity is then studied using these Hamiltonians on chains 
with periodic boundary conditions. The impurity is cou- 
pled to its first and second nearest neighbors (see Fig. [I]). 
For sake of simplicity, the couplings between two spins- ^ 
or between a spin-S impurity and one of its spin-i neigh- 
bors are taken identical, i.e. the same values of J, a, g 
and f2 are considered all along the chain independently 
of the nature of the spins. For instance, in the case of a 
spin-1 impurity, the Hamiltonian impurity part is explic- 
itly 

1, 



?C= J ]T(1 + <?(&, +6|))(SW 



(3) 



where the impurity spin operator is s and I (m) refers to 
its two nearest (next-nearest) neighbors. However, for a 
spin-0 impurity, the situation is simpler: all links between 
the impurity and its spin-i neighbors are taken as zero. 
In this case, when a — 0, the system corresponds to a 
chain with open boundary conditions (OBC). 




FIG. 1. Example of a periodic chain with L=5 spin-i sites 
and an additional impurity. Nearest neighbor (solid lines) and 
next nearest neighbor (dashed lines) links are plotted. Bold 
links correspond to the links between the impurity and its 
spin-i neighbors. 



B. Dynamical model 

Numerical technique: A reliable numerical treatment of 
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Hamiltonian (J]]) is a difficult task. Indeed, the phononic 
Hilbert space is infinite, even for a finite-jchain. Possi- 
ble approaches-such as QMC simulations^!, single-mode 
approximatioriE3_ or ED with a fixed maximum num- 
ber of phononsEa have been proposed. The results ob- 
tained in this paper are based on an ED calculation using 
the coherent states introduced by Fehrenbachcrta. Pre- 
liminary results have been reported in Ref. |l9| in the 
case of pure systems and Ref. ^l] for rings with im- 
purities in the absence of frustration. On each site, 
a two state basis including the phononic vacuum |0)j 
and a coherent state \l)i is considered. The phononic 
part of the Hamiltonian (Jy) can be rewritten as 7i p h = 
J A + b\)Ai + fl J^i , where A is a constant and 

Ai an operator of eigenvalues and 1 (the constant terms 
subtracted to Si ■ Sj in and (||) have been introduced 
with this purpose). It is convenient to introduce the co- 
herent state 

\l)i =exp(-r, 2 /2)exp (77 &J)|0><, 

where 77 is a variational parameter. In fact, is an 
eigenstate of the phononic Hamiltonian for fixed Ai = 1 
provided that rj = 770 = —XJ/Q. It is shown in Fig. ||(a) 
that T70 is the best value of 77 at small coupling. Fur- 
thermore it is a reasonable choice for arbitrary coupling. 
Therefore we shall assume 77 = r)o in the remainder of 
the paper. Note that the method could be improved by 
choosing explicitly the optimum 77 that minimizes the en- 
ergy for each set of parameters. The comparison between 
this coherent state approach (using 770) and a trunca- 
tion of the phononic Hilbert space retaining two or eight 
phonon states at each site is illustrated in Fig. ^(b) for 
the GS energy and the mean phonon number per site. 
The results are much more accurate using coherent states 
than when only keeping two phonon states even though 
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FIG. 2. (a) Relative GS energy of a 4-site chain as a func- 
tion of the variational parameter r\ for a = 0.5, g = 0.5, O = J 
(A) and a = 0, g = 1., O = J ((}). (b) GS energy (full line) 
and mean phonon number (dashed line) per site for a 4-site 
chain using a two (o), eight (•) phonon states truncation, and 
using the coherent state approximation with r\ — rjo (□) as 
a function of the magneto-elastic coupling g for a = and 

n = j. 



both Hilbert spaces have the same sizes. Note that al- 
though the GS energy is accurately obtained, the co- 
herent state basis seems to slightly overestimate the 
mean phonon number. However, even at large couplings, 
the phonon dynamics is qualitatively preserved. Con- 
sequently our non-perturbative approach takes properly 
into account the effects of the phonons and the results 
can be checked by comparing to a full phonon calcula- 
tion on small clusters (typically with L = 4 sites) as 
in Fig. ^(b). This variational treatment of the dynami- 
cal phonons hence enables us to handle clusters with up 
to L = 16 sites. Special care is needed for the treat- 
ment of the bonds connected to the spin-1 impurity. In- 
deed, special variational parameters (still corresponding 
to 770) have to be chosen for these bonds because the 
corresponding Ai operators for the phononic part of the 
Hamiltonian differ from those of the rest of the chain. 

Results: Let us summarize the main properties of this 
model. It has been shown that the GS of this model 
undergoes at T — a spontaneous symmetry breaking 
toward a dimerized gapped phase in a large region of 
parameter spaceE3. The elementary excitations are char- 
acterized as topological solitons, i. e. unpaired spins sep- 
arating the two different dimer patterns. Furthermore 
solitons (s) and antisolitons (s) do not bind. More explic- 
itly, one can define the ss binding energy as the p — > 00 
extrapolation of 



E S J(L = 2p) = [Eo(2p, 1) - E Q (2p, 0)] - 2e s , 



(4) 



where E (l, S z ) is the GS energy of a Z-site periodic chain 
in the S z sector. If I is even and S z = 0, E (l,S z ) 
corresponds in fact to the "vacuum" energy, the energy 
of the pure chain without any topological defect. e s 
is the soliton (antisoliton) minimum energy, defined as 
e s = lim p ^ oo [E (2p+l, ±)-£ *(2p+l)], where £ *(2p+l) 
is the interpolation between Eo(2p, 0) and Eg(2p +2,0). 
Thus, e s results as the energy difference between an odd 
chain containing a soliton and an (hypothetical) odd 
chain without soliton. Similarly, Eb{ss,L) is the energy 
difference between a chain with a ss pair and two iso- 
lated solitons. Previous studies provide strong evidences 
that the ss binding energy vanishesEj, i.e. Eg = in 
the thermodynamic limit in the absence of any explicit 
dimerization. Similar forms of binding energy will be 
discussed in the context of impurity-soliton binding. 

We first consider the case of a non-magnetic impu- 
rity. Typical patterns are shown in Fig. |§|(a) for parame- 
ters corresponding to CuGeOa. For even chains, the dis- 
tortion pattern 6i — g{b\ +b i ) rapidly oscillates with a 
two site period. The equilibrium pattern corresponds to 
strong bonds {i.e. Si positive) next to the impurity. The 
amplitude \Si \ varies along the chain with increasing mag- 
nitude for decreasing distance to the impurity. In other 
words, the impurity being located on site 0, the strongest 
bonds are the ones connecting sites 1 and 2 and sites L-l 
and L. In addition, the average magnitude of the modula- 
tion increases with a stronger magneto-elastic coupling or 
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weaker phonon frequency. The z-component of the spin 
(>S?) is zero because the chain is in a dimerized singlet 
state. For odd chains (Fig. ||(b)), the presence of a soli- 
ton at the center of the chain is signaled by a point where 
Si vanishes and (Sf) is maximum. This corresponds to 
the characteristic schematic picture of the soliton as an 
unpaired spin separating two dimer patterns. Similarly 
to the even chain case, strong bonds also form at the 
chain edges. The finite frustration that we have consid- 
ered here does not seem to change the main features of 
the pattern (see Ref. |2l] for comparison). The fact that 
the soliton forms at the center of the chain suggests that, 
contrary to the conclusions of Ref. solitons do not 
bind to non-magnetic impurities in pure ID spin-lattice 
models. 
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FIG. 3. Distortion Si (•) and z-component of the spin (S*) 
(o) as a function of the site i on (a) a 12-site and (b) a 13-site 
chain for a = 0.36, g = 0.8 and Q = J with OBC using ED. 

In order to be more quantitative, we define a soliton- 
impurity binding energy similarly to (||), 



E^iL = 2p + 1) = [E' (2p + 1, -) - E*(2p + 1)] 



e/, 



(5) 



where E' (l,S z ) is the GS energy in the S z sector of a 
chain with I spin-i sites and with an extra impurity, ej 
is an impurity minimum energy, ej = linip^oo [Eq(2p, 0) — 
Eo(2p, 0)], which indeed corresponds to the energy differ- 
ence between a chain with and without an impurity. 

The results shown in Fig. ^ (a) reveal no binding for 
any finite chain. Note that finite size effects are small 
for parameters in the vicinity of the MG point because 
of the small magnetic correlation length present in this 
case. 

The same calculations can be performed in the case of a 
spin-1 impurity. The different scaling behaviors in Fig. |] 
(pr for the spin-0 and j- for the spin-1 impurity) suggest 
that the physics involved is also different. Although a -jji 
behavior may also occur for larger rings with a spin-1 im- 
purity, our data are nevertheless compatible with a small 
binding between the spin-1 impurity and the soliton. 
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FIG. 4. Soliton-impurity binding energy using ED for 
a = 0.5, g = 0.4, Q = 0.3J (■) and a = 0, g = 0.45, Q. = 0.3 J 
(•) (a) as a function of the inverse square chain length in 
the case of a spin-0 impurity; (b) as a function of the inverse 
chain length j- in the case of a spin-1 impurity. 



C. Adiabatic model 

Numerical technique: As seen above, a full treatment 
of the lattice dynamics is tedious and restricted to small 
clusters. On the other hand, when f2 — > 0, an adiabatic 
approximation of the lattice is justified and enables us to 
handle larger clusters both in ED or QMC calculations. 
The main technical issue is to find the set of distortions 



{Si} which minimizes the total energy, 
leads to 



0, which 



(6) 



where (■ • • ) is the mean value in the GS for an ED cal- 
culation or some thermal average for a QMC calculation. 
Eq. (||) must be solved with the constraint of zero total 
distortion, 



5> = o. 



(7) 



This problem is clearly self-consistent since in the equi- 
librium condition (|^) the distortion pattern is present in 
both terms: implicitly in the first one because the GS 
depends on it, and explicitly in the second one. Conse- 
quently this probleiHnCan be treated using the following 
iterative procedure^ First, a randomly chosen set {<5°} 
satisfying (0) is taken. Applying the equilibrium condi- 
tion (p), where the mean values are considered using the 
adiabatic Hamiltonian with the set {<5°}, gives another 
set {5^} of distortions. After subtracting its mean value, 
a new set {S}} satisfying (0) is then obtained and the 
procedure is repeated until convergence. 

This scheme can be applied using ED or QMC tech- 
niques. Note however that in this last case, although 
larger lattice sizes can be considered, one is restricted to 
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the non-frustrated a = case to avoid the minus sign 
problem. As expected, for a (pure) periodic chain, the 
static uniform dimerization pattern Si = ( — l) l S is ob- 
tained through this procedure. 

Results: In this part, we shall study the introduction 
of a spin-0 or a spin-1 impurity in a chain described by 
the adiabatic Hamiltonian. It is here crucial to enable 
the magnitude of the modulation \Si\ to vary along the 
chain. A QMC self-consistent method is used so that we 
shall restrict ourselves to the non-frustrated a = case in 
order to avoid sign problems. A standard world-line al- 
gorithm was implemented. Most of the calculations were 
performed at T = 0.05J which we have shown in pre- 
vious studieSpis-Low enough to reflect essentially the GS 
properties H'tSEJ Besides, since we are mostly interested 
in GS properties we have carried on the simulations in 
the subspace of minimum z-component of the total spin. 
The number of slices in the Trotter direction was taken 
equal to M = 120. 

Patterns of dimerization are shown in Figs. ||and|| for 
a chain with a spin-0 impurity and are similar to those 
observed using the dynamical approach. For even chains, 
the main features are identical (Fig. |^). The distor- 
tion rapidly oscillates with a two-site period with strong 
bonds on the edges, where the impurity lies. Its am- 
plitude is slightly varying along the chain, the strongest 
distortion magnitude being observed at the edges. The 
magnitude of the distortion decreases when K increases 
as seen in Fig. ^. The z-component of the spin (Sf) van- 
ishes in this case because the chain is in the singlet dimer 
SP state. 
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FIG. 5. Distortion patterns Si on a 80 site chain as a func- 
tion of the position i for a — 0, K — J (•, the dotted line 
between the computed points is not shown for clarity) and 
K = 3J (o) using QMC (T = 0.05J). 

For odd chains, a single solitonic excitation is located 
away from the edges as shown in Fig. ^|. Its presence is 
inferred, as previously, from a local vanishing Si and a 
maximal (Sf). In this case we consider the local spin 
susceptibility, 



instead of (S?). Indeed, when using a QMC algorithm 
with an odd number of sites, a fictitious site has to be 
added and the total spin of the chain can fluctuate be- 
tween +i and — i through spin flips with this additional 
spin. However, the spin susceptibility Xi deals with a 
z-component of the lo cal spin with respect to the chain 
global spin orientatior£3ila. Note that, however, this ex- 
citation is topological in nature in the sense that Xi i n - 
tegrated in space on a finite region gives exactly a Curie 
law of a spin-i. The soliton width increases as the spin- 
phonon coupling decreases and becomes a sinusoidal dis- 
tortion in the weak coupling limit .til In fact, the inverse 
soliton width, the spin gap, and the amplitude of the 
dimerization are three different features of the SP transi- 
tion and all of them increase with the spin-phonon cou- 
pling. 
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FIG. 6. Distortion Si (a) and magnetic susceptibility \i 
(b) patterns as a function of the position i on a 79 site chain 
for a = 0, K = J (o, □ correspond to two different runs) 
and K = 2.5J (A) using QMC (T = 0.05J). Lines joining 
computed points are not plotted for clarity. 

It is interesting to notice that different QMC runs (o 
and □ in Fig. ||) lead to the very same solitonic patterns 
centered at different sites in a large area around the mid- 
dle of the chain. This is another evidence that there is 
no binding between non-magnetic impurities and solitons 
in this strictly one-dimensional model. The fact that the 
solitons are at least a distance away of the edges equal 
to half the soliton width indicates a likely short range 
repulsion between the soliton and the impurity. 

Next, we turn to the case of a spin-1 impurity. Typ- 
ical distortion and susceptibility patterns are shown in 
Fig. |?] for an odd length chain. In first approximation, 
the spin-1 inrpnrity and its two spin-i neighbors behave 
like a singlet,c3 i.e. like a spin-0 impurity as indicated by 
very strong bonds connected to the spin-1 impurity. This 
three-spin object is weakly connected to the rest of the 
system as indicated by particularly weak bonds (Si < 0) 
between sites 1 and 2, and L-l and L (the impurity is 
located at site 0)Jl3 Furthermore, strong bonds can even 
be observed close to this composite object (see arrows in 
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Fig. 0) as it is the case for a real S=0 impurity. Physi- 
cally, the three-spin cluster spends most of the time in a 
singlet state configuration. One observes that the soliton 
always stays quite close to the impurity, independently 
of the QMC run, suggesting a small residual impurity- 
soliton binding, consistent with the previous dynamical 
calculation. 
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FIG. 7. Local distortion Si (•) and z-component of the local 
spin (S^) (o) as a function of the position i on a 79 site chain 
for a = and if = 1.5 J using QMC (T = 0.05 J). A spin-1 
impurity is located at site at the middle of the figure and 
periodic boundary conditions are used. Two arrows indicate 
particularly strong bonds (Si > 0) next to the spin-0 object 
formed by the spin-1 and its two neighbors. 

III. EXPLICITLY DIMERIZED CHAINS 

The Hamiltonians considered in Sec. [njare purely one- 
dimensional. However, in SP compounds, the phonons 
have a three-dimensional character. One of the sim- 
plest models proposed in the literature to take this 
three-dimensional behavior into account is the dimerized 
Heisenberg chain, 

Hmf = .7 + ■ S i+1 + aS t ■ S i+1 ], (9) 

where S is an explicit dimerization. Contrary to the mod- 
els considered in Sec. the distortion pattern is here 
static and uniform. It is interesting to note that this 
Hamiltonian (|9|) can in fact be inferred from the dynam- 
ical Hamiltonian (Fni— when only the dominant 7r-mode 
phonon is consideredcll. 

The dimerized Heisenberg chain has been widely stud- 
ied in the literature. Its main features are the following: 
the elementary excitations can be interpreted as mag- 
netic solitons and antisolitons (or spinons), namely isn- 
lated spins- i separating two dimerization patternsQ'EH. 
As in the models considered in the previous section, for 
even site chains, solitons and antisolitons do not exist in 
this model as independent particles but appear in pairs. 
However, in this difference to those models. 



they are confined by a linear., potential proportional to 
8(^= 0) in the weak 6 limitCJ'OE3. It has also been shown 
that the excitation spectrum is constituted of a ladder 
of soliton-antisolit on .. bou nd states lying below a two- 
magnon continuumc!Jl23EZl. The differences between the 
models studied in the previous section and the present 
one can be schematically visualized in the large dimer- 
ization limit as in Fig. |^. While in the fully quantum dy- 
namical model both the distortion pattern {5i} and mag- 
netic dimers are two-fold degenerate, in the chain with 
fixed dimerization only one of the distortion patterns is 
chosen and a magnetic state with ss pair separating two 
magnetic dimers configuration would leave a region with 
a "wrong" pattern in between the soliton and anti-soliton 
thus leading to a confining potential. On the other hand, 
this confining is absent if Hamiltonian (|^) is considered 
since in this case the lattice would relax following the 
magnetic order. 
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FIG. 8. (a) The two-fold degenerate GS of Hamiltonian 
([j]). The solid (dashed) lines indicate strong (weak) bonds 
and dimer singlets, (b) The fixed dimerization pattern of 
Hmf (above) and a magnetic configuration of singlet dimers 
with a ss pair (below), (c) A magneto-elastic state without 
(above) and with a ss pair (below) of Hamiltonian (^) . 

The case of an odd chain with open boundary con- 
ditions, i.e. with a spin-0 impurity cutting the chain, 
has been extensively studied in Ref. [l^. It has been 
shown that, due to the linear confining potential, an iso- 
lated) s o li ton binds to the spin-0 impurity next to a weak 
boncfcjOEa. In this section, we generalize the approach 
of Ref. |l5] to investigate the possibility of several bound 
states and to study the low energy spectrum. Following 
Ref. [56], we assume that the dimerization is stabilized by 
a strong frustration a ~ 0.5 while the "confining force" 
introduced by S remains small. The excitation spectrum 
with respect to the "vacuum" energy i.e. the GS energy 
of a periodic even chain (without any defect) interpolated 
to odd chains, is shown in Fig. |^. The energies could also 
be equivalently measured with respect to the beginning 
of the continuum but finite size effects affect its posi- 
tion making results more dependent of the lattice size. 
For a = 0.5 and S = 0.05, four soliton-impurity bound 
states below the continuum can be clearly identified. The 
energy of the soliton-impurity bound state can increase 
up to an upper bound, above which it becomes energeti- 
cally more favorable to create a soliton-antisoliton bound 
state. This state, which is a combination of a soliton- 



G 



impurity and soliton-antisoliton bound states, has a spin 
| and is the lowest state of the continuum. We have 
checked that its energy £7 (§) = (0.7949 ± 1CT 4 ) J, is in 
the thermodynamic limit the sum of the lowest impurity- 
soliton bound state E {\) = (0.2748 ± 1(T 4 ) J plus the 
lowest energy of a ss bound state, i.e. the spin gap A 01 = 
(0.5200 ± 10~ 4 ) J, where A 01 = E (2p, 1) - E Q {2 Pl 0) in 
the notation of Eq. (|4|). The above energies are infi- 
nite size extrapolations obtained with ED and DMRG 
for a = 0.5 and 5 = 0.05. 




1/1/ 

FIG. 9. Lowest lying S z = \ (o) and S z = § (X) excita- 
tions for the q = 0.5 and S = 0.05 Heisenberg chain as a func- 
tion of the square inverse of the chain length L (up to L=25) 
obtained by ED. The energy reference is the GS of the even 
periodic chain interpolated to odd sizes. Dashed lines repre- 
sent a {-£z, T3-) fit and the extrapolations to infinite sizes are 



indicated. DMRG extrapolations are also plotted (♦) 
dot-dashed line indicates the onset of the continuum. 



The 



The wave functions (Sf) obtained in a DMRG calcula- 
tion shown in Fig. mi clearly support this scenario. The 
four lowest states (a^b,c,d) show a soliton bound to the 
weak bond edge, and the soliton moves further away from 
the impurity when its energy increases. In contrast, in 
the lowest spin-| state (Fig. [H](e)), a ss bound pair can 
be clearly identified in addition to the solitonic GS seen in 
Fig. |l^(a). The number of soliton- impurity bound states 



below the continuum seems to increase with 



s was 



originally proposed for the number of ss pairscJ 



IV. ELASTIC INTERCHAIN COUPLING 

The three-dimensional character of phonons has been 
considered in Sec. [II using the dimerized Heisenberg 
chain. However this can be achieved in a more appro- 
priate way by adding to the pure ID models of Sec. || 
an interchain coupling. Such an approach is particularly 
necessary for materials with large anisotropy in the elas- 
tic properties. Our motivation here is to improve the 
description of Sec. Ill in order to enable the lattice to 



locally relax following the magnetic order. As seen pre- 
viously in the case of isolated chains, such effects are 
crucial in the presence of inhomogeneities introduced by 



local magnetic excitations (solitons) or impurities and 
lead to important qualitative differences with the dimer- 
ized Heisenberg chain (||). 
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FIG. 10. {Sf ) as the function of the site i for the four lowest 
spin-^ (a,b,c,d) and the lowest spin-| state (e) for a — 0.5, 
S — 0.05 on a 51-site chain using a DMRG algorithm. The 
bond between the first two sites is "weak" . In (e) , a magnifi- 
cation by a factor 6 has been applied to the points left to the 
dashed line. The 7 first sites of the lowest bound state (a) 
(A) and a soliton-antisoliton pair on a 44 site lattice (•) are 
also shown for comparison. 

In a quasi-lD SP compound, a given chain cut at its 
ends by spin-0 impurities ( "impurity chain" ) is immersed 
in the bulk. The neighboring chains are in the SP phase, 
i.e. they have a uniform dimerization, and produce on 
this chain a q = ir po tent ial through an elastic interchain 
coupling of the formt§E3, 



where Si is the distortion of the impurity chain, {5f xt } are 
associated to the neighboring chains and K± is the in- 
terchain elastic coupling constant. Note that using this 
model, two neighboring chains can be in phase or out 
of phase depending on the sign of Kj_. The neighbor- 
ing chains can be treated in the mean-field approxima- 
tion while fully retaining the dynamics of the considered 
chain: 5° xt = (— l) J+1 <5o- Consequently, one obtains, 



i+l j 



where K' = K±Sq. 

As schematically shown in Fig. ll|(a), in a strong exter- 
nal potential, the situation is similar to the one of Sec. Ill 
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the soliton is bound to the spin-0 impurity at the weak 
bond edge. However, in a weak external potential, results 
of Sec. H suggest that solitons can delocalize away from 
the spin-0 impurity edge as depicted in Fig. |ll|(b) leaving 
strong bonds at the edges. Our aim here is to study how, 
when switching on an external potential and increasing 
it, the soliton binds to the impurity. The results shown 
in this part are obtained using an adiabatic treatment of 
the lattice and a QMC si mula tion of the spin degrees of 
freedom (similar to Sec. 1IC) for non-frustrated chains 
with spin-0 impurities. Consequently, the lattice can be 
seen as open chains with two unequivalent non-magnetic 
impurities at its edges. 



(a) 



(b) 



— o t — — — o — — t 



FIG. 11. Introduction of a non-magnetic impurity in a 
three-dimensional spin-Peierls compound. Schematic pattern 
pictures are drawn for the case of K' > 0. (a) corresponds to 
a strong external potential and (b) to a weak external poten- 
tial. Full (dashed) lines symbolize strong (weak) bonds. 

We first consider the case of an odd size chain. For van- 
ishing K', the soliton is free and strong bonds (Si > 0) 
occur at the two chain ends. A non-zero K' creates a lin- 
ear confining potential between the soliton and the impu- 
rity where the potential tends to form a weak bond. This 
confining potential then equilibrates with the short range 
impurity repulsion. Consequently, the soliton moves to- 

0.2 
0.1 

-0.1 
-0.2 
0.1 

to - o.o 

-0.1 
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to - o.o 

-0.1 
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-0.1 - 
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FIG. 12. Lattice distortions Si as a function of the position 
i on a 79 site chain for a = 0, K — 2 J with an interchain 
coupling K' obtained by QMC simulation (T = 0.05 J). The 
value of K' is indicated on the graph. 



wards the impurity when the interchain coupling K' in- 
creases, as can be seen in Fig. [l2|. Moreover, the magni- 
tude of the dimerization increases with K' as expected. 
Note that the distortion pattern obtained in Fig. |l2|(d), 
i.e. stronger bond on one edge and zero distortion at 
the other edge, is strikingly different from those of Rcf. 
p0| . These two different behaviors (obtained for i ~ 1 
and i ~ L) are in fact expected on each side of any given 
impurity in this system. 

For even size lattices, two cases have to be examined. If 
the free distortion pattern and the external potential are 
in phase (i.e. if K'S^Sf* > 0, where <5° = 5 l (K' = 0) ), 
the external potential amplifies the free distortions seen 
in Fig. [| with strong bonds on the edges. However, if 
they are out of phase, the external potential will create 
two topological soliton defects. These excitations move 
to the chain edges and become more localized when the 
external potential increases; .similarly to what happens 
for a static dimerized chainE3C3. The patterns shown in 
Fig. [l3 support this scenario: a spin ±~ is located near 
each edge. Note that at low temperature, this excess of 
spin-i at the chain's ends oscillates between positive and 
negative value for different QMC runs, corresponding to 
the S z — > — S z symmetry. Consequently the role of the 
interchain elastic coupling is crucial in order to stabilize 
impurity-soliton bound states. 
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-0.2 
A 0.1 
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FIG. 13. Lattice distortions 5i (•) and z-component of spin 
(Sf) (o) as a function of the position i on a 80 site chain for 
a — 0, K = 2J with an interchain coupling K' using QMC 
(T = 0.05 J). The value of K' is indicated on the graph. 

The formation of the ss topological defect in the open 
chains with an even number of sites L (as described 
above) is particularly subtle. In fact, because of the 
finiteness of the energy cost A 01 (or 2e s if the solitons 
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are far away) associated with the formation of the de- 
fect, one can, on general grounds, deduce the existence 
of a critical value K' c of K' in such a phenomenon. A 
comparison between the energy cost oc 2e s ~ A 01 and 
the transversal elastic energy gain oc K'(L — 2T), where 
T is half the soliton width, leads to K' c ~ A 01 /(L - 2T). 
Since A 01 increases with decreasing L ) and T = T(K) 
is roughly independent of L (see Ref. [58|), then we ex- 
pect that K' c decreases for increasing chain length. For 
K = 2 we have obtained K' c = 0.050,0.032 and 0.005 
for L = 20, 40 and 80 respectively, thus confirming this 
prediction. 

For real SP compounds doped with non-magnetic im- 
purities, chains are cut at random places. By simple in- 
spection of all possible configurations, the previous study 
shows unambiguously that exactly one spin-^ becomes 
bound to each impurity due to the interchain elastic cou- 
pling responsible for the three-dimensional character of 
the lattice dynamics provided that K' > K' c . How- 
ever, in chains short enough that K' C (L) > K', non- 
magnetic impurities at the ends do not carry magnetic 
moments. In other words, let us assume that one, by 
thought, moves closer to each other two impurities with 
respectively a soliton and anti-soliton bound to them. 
Our previous analysis shows that, below a critical sepa- 
ration L c ~ 2r+A 01 /if ' between the two impurities, this 
magnetic state becomes unstable towards a new config- 
uration where the soliton and the antisoliton annihilate 
each other and leave behind a dimerization pattern in an- 
tiphase with the 3D order. Notice that such open chains 
are most likely to appear in the system after doping since 
the probability distribution of lengths is P(L) — c(l— c) L , 
where c is the impurity doping. Furthermore, the char- 
acteristic distance L c can become fairly large for small 
interchain elastic coupling K 1 . 

Finally, in Fig. [li] we show the staggered spin- 
spin correlation function C^ s \\i — j\) = (1/£)X^ j (Si ■ 

£/>(- 1 ) l< " il (normalized such that C^(0) = 1) for 
L = 80, K = 2 and K' = and 0.02 obtained by 
QMC. A clear enhancement of the AF correlations can 
be seen when a soliton-impurity bound state is present 
(K 1 = 0.02) in contrast to the case of an almost uni- 
form pattern with strong end bonds shown in Fig. |B| 
{K 1 = 0.0). This result is consistent with the enhance- 
ment of AF correlations close to an impurity for weak 
end bonds found in Ref. [l5| One expects that this in- 
crease of the AF correlations in the vicinity of impurities 
carrying an effective spin-1/2 magnetic moment can be 
observed in inelastic neutron experiments at low energy 
transfer and momentum transfer q ~ n. 

Our results are also consistent with the experimen- 
tal observatioruof Curie laws in the susceptibility of Zn- 
doped (spin-O)Ej and Ni-doped (spin-l)EJ CuGeOa ma- 
terials. Note however that, in the case of Zn-doping, due 
to the mechanism described in Sec. |y|, a significant frac- 
tion of the impurities may not carry an effective spin-1/2 
leading to a reduction of the magnitude of the Curie 



0.20 p^r 




FIG. 14. Staggered spin-spin correlations in the 80 site 
chain with K = 2 obtained by QMC (T = 0.050). Peri- 
odic boundary conditions with K' — and OBC with various 
values of the interchain coupling K' as indicated on the plot 
have been considered. 

term. In the case of Ni-doped CuGeOs (Ref. [59]), an 
interesting cross-over at low temperature was observed 
between two different Curie law behaviors. We interpret 
this low temperature behavior with the formation of an 
effective singlet at the Ni-site which can bind a spin-i 
soliton as described in Sec. ffl. 



V. CONCLUSION 

In this work, we have studied several spin-lattice 
models using various numerical techniques. Local non- 
uniform dimerization patterns have been obtained result- 
ing from inhomogeneities due to impurity doping. We 
have shown that the lattice responds strongly to the vari- 
ation of the spin density even in the adiabatic case. Bind- 
ing between elementary solitonic excitations and spin-S 
impurity (S=0 or 1) has been investigated. In strictly ID 
models, no binding between solitons and non-magnetic 
impurities occurs. However, a small binding between 
spin-1 impurities and solitons has been inferred. The 
role of the three-dimensionality of the lattice was also 
investigated by comparing a model with a fixed mean 
field-like dimerization to a model with an explicit cou- 
pling to an adiabatic lattice which can locally relax to 
follow the magnetic order. In the latter case, the inter- 
chain clastic coupling was shown to be responsible for the 
binding of solitons to spin-0 impurities. However, in con- 
trast to the dimerized Heisenberg model, we predict in 
this model the existence of S = impurities carrying no 
spin-1/2 solitonic excitation. These impurity sites must 
appear in pairs separated by sufficiently small chain seg- 
ments with an even number of sites and a dimerization 
pattern in antiphase with respect to the 3D dimerization 
order. 
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